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Traditional quantum error correction involves the redundant encoding of k quantum bits using 
n quantum bits to allow the detection and correction of any t bit error. The smallest general t = 1 
code requires n = 5 for k = 1. However, the dominant error process in a physical system is often 
well known, thus inviting the question: given a specific error model, can more efficient codes be 
devised? We demonstrate new codes which correct just amplitude damping errors which allow, for 
example, a t = I, fc = I code using effectively n = 4.6. Our scheme is based on using bosonic states 
of photons in a finite number of optical modes. We present necessary and sufficient conditions for 
the codes, and describe construction algorithms, physical implementation, and performance bounds. 



I. INTRODUCTION 

Information is classically measured in units of bits, 
which are deterministic two-state systems that are said 
to exist either as logical zero, Oz, or logical one 11- How- 
ever, such a representation is only an approximation of 
reality, which is described at small size scales by the laws 
of quantum physics. Quantum-mechanical two-state sys- 
tems ( "qubits" ) differ from the classical bit in that they 
may exist in a superposition of the two states, for ex- 
ample as = a\Qi) + ^IIl); where a and h are arbi- 
trary complex coefhcients which satisfy |ap -|- |&p = 1. 
Two continuous real parameters are needed to describe 
the state of a qubit, and in this sense, more information 
is somehow carried in it than by a classical bit. Fur- 
thermore, qubits may not be cloned [0, and even more 
importantly, they may exist in entangled states where, 
for example, two qubits only carry one qubit of quantum 
information. 

Unfortunately, quantum information is (partially) lost 
whenever a quantum system is observed (whether delib- 
erately or inadvertently) . This decoherence process plays 
a role analogous to noise in a classical communication 
channel. A major advance in quantum information the- 
ory has been the discovery that quantum information can 
be redundantly encoded in such a manner that it may 
be efficiently transmitted with arbitrarily high fidelity 
through a decohering quantum channel. The quantum 
error correction codes [^|j3j which make this possible are 
analogous to classical codes for binary memoryless chan- 
nels. Corresponding codes for classical linear codes and 
Reed-MuUer codes have been found |^,||. 

These quantum coding schemes are based on a model 
for the decoherence of qubits, in which three kinds of er- 
rors can occur: bit fiips (|0) ^ |1)), phase flips (|1) ^ 
— 11)), and both simultaneously. This model is general; 



it describes all possible decoherence mechanisms for a 
qubit. However, in a given physical system, the dominant 
decoherence process is of a specific nature which may ad- 
mit a simpler description. For example, in phase damp- 
ing, no bit flips occur. The question therefore arises: 
given a particular decoherence process, what is the opti- 
mal quantum error correction scheme? 

We do not yet know how to handle this general prob- 
lem. However in this paper, we report on progress to- 
wards a solution, by demonstrating a new class of quan- 
tum error correction codes which correct only one partic- 
ular decoherence process known as amplitude damping. 
Our approach is similar in philosophy to that of 1^ , but 
in contrast to that and other previous schemes, instead of 
qubits, which live in a two-dimensional Hilbert space, we 
utilize bosonic systems which occupy the Hilbert space 
|0) • • • |iV). We are unaware of any classical analog to our 
codes. We present possible physical implementations of 
our scheme, and conclude with a comparison with exist- 
ing binary codes. 

II. AMPLITUDE DAMPING MODEL 

Noise is a fundamental process which accompanies the 
dynamics of any open system. Traditionally, the dynam- 
ics of an open quantum system are described by a "mas- 
ter equation" The system a, described initially by 
the density matrix pao, couples to a bath h through an 
interaction Hamiltonian Hi. Evolution generates an en- 
tangled state of the total system p, and we trace over the 
bath state at the end. Taking the system to be a simple 
harmonic oscillator, and using the simplest bilinear in- 
teraction in which photons are exchanged back and forth 
with the bath, we find in the Born and Markov approxi- 
mations an equation of motion for the density matrix 
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2apa') . 



(2.1) 



where A is a coupling constant between the system and 
the environment. We have set (6^6) = in (2.1) to reflect 
a reservoir at temperature kT much smaller than the sys- 
tem's energy scale huj. This process describes the gradual 
loss of energy from the system to the environment, and 
is known as amplitude damping 

Mathematically, the evolution of the density matrix 
from t to t' in a particular process may be described 
as a linear transformation from one density matrix p to 
another, p' . This may be expressed (in the "positive op- 
erator sum representation" ) as 



k 



AkpAl , 



(2.2) 



where A/j are positive linear operators, sometimes re- 
ferred to as "Krauss effects pl," which are related to the 



Lindblad operators appearing in Eq.(2.1). For amplitude 
damping, we find that 



(2.3) 



where 7 = 1 — cos^ x is the probability of loosing a sin- 
gle photon from the system between the initial and final 
times. After tracing over the final environmental state, 
we find the operators in the Hilbert space of the system 
alone: 



A, 




k) {n\ 



(2.4) 



If the initial state is pure, it may be written as p = 
The final state p' may be elegantly described as 
an explicit mixture of pure states given by 



N 



(2.5) 



fc=0 



where N is the maximum occupation number of a single 
bosonic mode. Here, the "©" symbol represents a tensor 
sum of states, and [-0') is a convenient shorthand used to 
denote a mixed state, as distinguished from a pure state 
\ih). In other words, 



N 



(2.6) 



The mixed state is a tensor sum of + 1 (un- 
normalized) pure states which describe the + 1 pos- 
sible final states of the system; one may interpret these 
as non-interfering "alternative histories p^." The nor- 
malization of each pure state gives its probability of oc- 
currence. In general, k describes the number of photons 
lost to the environment. It is important that even when 



no photons are lost to the environment, then the state of 
the system is changed. 

So far, we have described the effect of amplitude damp- 
ing on a single mode system. Consider now a system with 
m modes, and let us use A^j to denote the action of the 
effect Ak on the j*'* mode of a state, j £ [l,ni]. After 
amplitude damping, the initial pure state 



becomes the mixed state 



[i>out) = 



N 

0^feo|?^o) 

k=0 



N 



Akn 



.k=0 



(2.7) 



(2., 



where there are now (A^ + 1)™ possible final states. It is 
convenient to use the shorthand notation 



Az 



Akoo ■ 



(2.9) 



where kj is the j digit of the number k written in base 
N-l-1, so that we may rewrite Eq.(2.8) as 



[ipou 



fc=0 



(2.10) 



Note that identical states in a tensor sum can be com- 
bined using the rule 



(2.11) 



since an overall phase doesn't matter (assuming no en- 
tanglement with other systems). 

As an example, amplitude damping of the state 



|V'^n) =a|01) +&|10) 



gives, using 



^o = |o)(o| + vT^|i)(i| 

^i-V7|0)(l|, 



(2.12) 



(2.13) 
(2.14) 



the output state 

i^out) ^ Aaol^Pin) ® Aoil^in) ® Aiol^P.n) ® Aiil^in) (2.15) 

(2.16) 



This result can be understood intuitively: the original 
state only contains a single photon, and thus, whenever it 
is lost, the final state must be the vacuum. This example 
indicates that the state of Eq. ( 2.12| ) is useful for detection 
of a single photon loss. However, no useful information 
about a and b can be extracted from the vacuum state, 
and so it is not useful for error correction. 
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III. EXAMPLES 



with the most probable state: 



Let us motivate the remainder of this paper by consid- 
ering the following example: We encode the logical zero 
and one states of a single qubit as 



|0l> 



1 40) + 1 04) 
V2 



II; 



|22) 



such that the initial state is the arbitrary qubit 
IV-m) =a|0L)+5|lL). 



(3.1) 



(3.2) 



The possible outcomes after amplitude damping may be 
written as 



(3.3) 



where we shall express fc as a base 5 numeral, and |(/)j,) is 
an unnormalized pure state (the norm of which gives its 
probability for occuring in the mixture). For small loss 
probability 7, the most likely final state will be 



I0OO) = (l-7)'l^m) 



(3.4) 



corresponding to no quanta being lost to the bath. The 
next most likely states result from the loss of a single 
photon: 



27(1-7)3/2 L|03) +6|21) 



|0io) = V27(l - 7)'^' a|30)+6|12) 



(3.5) 
(3.6) 



States resulting from the loss of more than one quan- 
tum occur with probabilities of order 7^. Therefore, we 
limit our correction scheme to errors losing at most one 
quantum. Each such error Ei takes |0l) and to 
states IOl)^ and \lL)i respectively. The key is that |0l), 
|1l)j \^L)i and {Ihji Vi are mutually orthogonal, and 
so are |0oo)j I0oi)j and \4>w)- In principle, a ("quan- 
tum non-demolition") measurement scheme can detect 
all error syndromes. Furthermore, for each i, the norms 
of IOl)^ and |lL)i are equal. After detecting an error 
syndrome, one can apply an appropriate unitary trans- 
formation converting \QL)i and |lL)i, to |0l) and |1l) 
respectively. This makes possible the correction: 



a|Oz 



-b\lL)^^ a[a\QL) +h\lL)] 



(3.7) 



where a is independent of a, h. Note, this is done with- 
out any information about a, 6, and without diminishing 
the amplitude of the erroneous state. For this particu- 



lar code, the output state has fidelity (see Eq.(|7.3|) |ll[ ) 
J- = 1 — 67^ with respect to the input. 
As a comparison, consider the code: 



|0z 



111) 



11/ 



|22), 



(3.8) 



|0oo) = a(l - 7)|11> + - 7)'|22) . (3.9) 

No unitary transformation will bring it back to 

a|ll)+6|22). (3.10) 

unless a, h is predetermined (a non-unitary transforma- 
tion can revert the change, but it will reduce the fidelity 
of the correction process to 1-0(7)). 

In the remainder of the paper, we shall describe the 
criteria for a scheme in which k qubits may be encoded 
so that loss up to t quanta may be corrected. For small 
t, a scheme will be exhibited. 



IV. CODE CRITERIA 

Quantum error correction is just the reversing of some 
effect due to decoherence. General criteria for this to 
be possible have been given in the literature . In this 
particular case, we may express the required conditions in 
the following manner. Let {|co) • • • |c;) • • • jc;^)} be /q + 1 
codewords which encode orthogonal logical states within 
the m mode Hilbert space with maximum total photon 
number N , and define JC{t) as the set of all m-digit base 
N -\- 1 numbers whose digits sum to t (corresponding to 
t errors). The logical states must satisfy 



{ci,\A\Ay\ci,)^0 



{ci\A\A^M) 



9k 



for h^hoT k' (4.1) 
Vi (4.2) 



for all fc, fc' e Us<t ^(s)- Here, gj^ is some constant which 
depends only on fc. This is an extension of the non- 
degenerate quantum error correction code criteria given 
by ll^] to the case where none of the Ak operators are 
identity. The first condition requires that all erroneous 
states be orthogonal, and the second requires that the en- 
coded Hilbert space not be deformed. Here, we present 
an explicit statement of these two conditions as algebraic 
conditions on the code construction. 

We consider first a codeword \ci) which may be ex- 
pressed as an equally weighted sum of A^; energy eigen- 
states, \ni...nm} (we shall refer to these as quasi- 
classical states, "QCS" for short). When all the QCS 
are equally weighted, we call the code "balanced" . Oth- 
erwise, the code is referred to as "unbalanced". Each 
codeword can be represented by a matrix with m columns 
and Ni rows, each row being one of the QCS in the code- 
word. For instance, if 



1 



iq; = 



|"11 • ■ • ?^lm) + • • . + \nNil ■ ■ ■ nN,m.) 



(4.3) 
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then the corresponding matrix Mi is: 



nil ni2 

"-21 n22 

nNii nNi2 



nir, 

n2-n 
nNvn 



For t = 
Ao\nii ...riy 



errors, we have /C(0) 
.> = (l-7) 



RSi/2 



mi 



(4.4) 



= {0}, and 
where the row 



sum RSi = X^i'Li "^ij- Criteria Eq.(4.2) requires the am- 



pUtudes of AqIq) be the same for ah |q), that is: 



1 ^' 



i=l 



be the same for all |q). A sufficient condition for this is 
the equality of all the RSi for all i and for all codewords 
\ci). Alternatively, one may say that the sum of any row 
from any JVli equals TV, the same total photon number in 
all the QCS. Denote the set of QCS with m modes and 
total photon number N as Q{N,m). It follows that if 
we construct all the codewords from states in Q{N,m), 
then the non-deformation constraint Eq.( |4.2| ) is satisfied 
for i = 0. Physically, this requirement stems from the 
fact that a state with higher number of quanta decays 
faster. To preserve the a posteriori probabilities of each 
codeword, we must encode them in a subspace in which 
the decay probabilities are equal for all of them. 

Similarly, for t = 1 errors, we have /C(l) — 
{0 • • • 1, • • • 10, . . . , 1 • • • 0}, and, for example. 



Ni 



Aq...i\ci) = Ao...i-^ V |nii . 

^=l 



(4.6) 



N, 



2_^\ rr \na ■ ■ -riim - 1} , (4.7) 



and 



(q|^J...i^o-i|q) 



Ni 



7(1-7)^-^^, 



N, 



i=l 



(4.8) 



The non-deformation criteria requires the above sum to 
be the same for all codewords. Equivalently, the column 
sum of the to*'' column of each codeword divided by Ni 
has to be independent of the codeword. Similar expres- 
sions for other A^, give rise to similar criteria for each 
column separately. 

We therefore have the following: 

Lemma 1.1: Let each codeword be expressed 
as an to column, Ni row matrix with el- 



ements n. 



■ij- 



If we choose codewords such 



that Y^i^ijl^i 



(Q|4A^,h) = <7fc,VAe/C(l) 
Proof: As above. □ 



Hj for all |ci), then 



This criteria corresponds to certain symmetry require- 
ment among the various codewords. For t = 2, 
/C(2) = {0-- •02,0- ••11,0 •••20...,20-- -0}. Working 
out Aj,\ci) for each k and applying the criteria for t = 
0, 1, 2, one arrives at the following: 

Lemma 1.2: Same setting as LeTOTOol.l. Let 
us choose codewords which satisfy the non- 
deformation criteria for t — \, and such that 
J^rntj.UijjNi = yj^j^ for aU |q), where 
h,j2 G [1, A^i], ji,j2 niay or may not be dis- 
tinct. Then {cl\A\AJ^,Cl) = Vfc e /C(2). 

Proof: See Appendix A. □ 

Generalization to arbitrary t is as follows: 

Theorem 1: Let each codeword be expressed 
as an to column, Ni row matrix with ele- 
ments . If we choose codewords such that 



%ij2,...j. indepen- 
.,j,) G [1,A^,]^ and 
Vfc G 



9h 



^2 ^ijl "ij2 ' ' ' ^ijt I -^l " 

dent of |ci), V/, ^{31,32, ■ ■ ■ 
Vs e [l,t], then {ci\AIA 

Proof: See Appendix A. □ 



The above theorem can be generalized to unbalanced 
codes in which codewords are unequally weighted su- 
perpositions of QCS. If the amplitudes of the QCS 
in \ci) are {^/JIi, ^/Jja, ■ ■ ■ , ^J^'Ni ) , we replace the sum 

„ • • • rii,,, i.e., we 



it 1^1 by 



replace the equal weights by the ^i's (the derivation 
of Theorem 1 in the unbalanced case is a straightforward 
generalization of the balanced case and we will skip the 
proof). 

As t increases, the non-deformation criteria becomes 
very restrictive. We have found unbalanced codes by nu- 
merical search correcting up to i < 4, as (Section VIII) 
which have no analogues in the balanced codes. On the 
other hand, for t < 2, we found simple construction algo- 
rithms for balanced codes with no apparent counterparts 
for the unbalanced codes. 

It should also be noted that the i = non-deformation 
criteria, that row sums (total number of excitations in 
each QCS) be equal for all rows and for all matrices 
(codewords), is not a necessary condition. An example 
is Shor's 9-bit code [|: 



|0l)=(|000) + |111))^^ 
|1l)=(|000)-|111))«^ 



(4.9) 
(4.10) 



The QCS have different number of I's, but Eq.(4.5) is 
equal for jO^) and |1l)- However, code criteria for t > 1 
will be extremely complicated when row sums are differ- 
ent, and treatment of such codes are outside the scope of 
the present discussion. 
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The other criteria, the orthogonality constraint 
Eq.(4.1) can be satisfied as foUows. Let \u) = \ui . . . Um) 
and \v) = \vi . . . be two states in Q{N, m). We define 



the distance between u and v as 



(4.11) 



Clearly, Q < V < N. Moreover, V{u,v) = V{v,u), 
'D{u, u) = 0, and 

V{U, V) + V{V, w)^^^\u,-V^\ + \V^ - W^\ (4.12) 

i 

i 

= V{u,w). (4.14) 



Thus I? is a metric on the discrete space Q{N,m). (For 
binary states, this is half of the Hamming distance). De- 
fine the distance between two codewords |ci) and |c2) to 
be the minimum of I?(wi,U2) with \u2) being QCS 
in |ci) and |c2) respectively. For codewords with non- 
negative amplitudes of the constituent QCS, two code- 
words are orthogonal iff their distance is non zero. We 
therefore have the following: 

Theorem 2: Let |ci) and |c2) be two code- 
words formed from states in Qi and Q2 re- 
spectively, where Qi,Q2 C Q{N,in) and 
2?(mi,M2) > t Vui £ Qi,U2 G Q2- Then 
{ci\AlAj^,\c2) =0,yk,~k' e[j,^,IC{s). 

Proof Let Ap\ci) ^ \di), A^|c2) = ^2), 
and let \vi), \v2) be QCS in \di), 1^2) respec- 
tively s.t. T>{di,d2) = 'D{vi,V2)- Let \ui), 
1^2) be the original QCS in |ci), |c2) before 
the error. Then, I?(7ii,wi) = 'D{u2,V2) < 
t/2, and V(vi,V2) +V{ui,v{) +V{u2,V2) > 
2?(mi,M2) > t. Hence, 'D{vi,V2) > and 
V{di,d2) > 0. Therefore, \di) and 1^2) are 
orthogonal states. □ 

In other words, by forming codewords using QCS which 
are sufficiently far apart, then the orthogonality condi- 
tions is easily satisfied. 



V. EXISTENCE OF CODES 

How large must N, m, and Ni be to satisfy both the 
non-deformation constraint, Eq.(4.2), and the orthogo- 
nality constraint, Eq.(4.1)? Wc now show that an unbal- 
anced code exists for arbitrarily large t if is allowed 
to be arbitrarily large, and give an upper bound for the 
required N. 

Let |co), |ci), |q), |c;J, helo + l codewords, 
each being an unequally-weighted superposition of Ni 
QCS in Q(n,m). For convenience, define 



V{n,m)=C{n + m-l,m-l)= + (5.1) 

n!(m — 1)! 



as the number of all possible partitions of the integer n, 
i.e., the number of xi such that xi + X2 + ■ ■ ■ -\- Xm = n 
ilii . Then if we choose N = nd such that 



Tin, m) > TVi + 7V2 + • • • + A^/o = Nqcs ■ 



(5.2) 



where Nqcs is the total number of QCS in the code- 
words, then by Theorem 2, all the QCS involved can be 
chosen to be distinct, and multiplication of the number 
states by d = t + 1 allows the orthogonality condition to 
be satisfied. 

On the other hand, the non-deformation condition in- 
volves satisfying a certain number of constraint equa- 
tions, given by the total number of possible errors times 
lo- The number of errors involving losing s photons from 
m modes is just the number of partitions of s into m 
parts, V{s,m). Take the QCS to be arbitrary, and solve 
the non-deformation constraint equations (of Theorem 1, 
generalized to include unbalanced codes) as linear equa- 
tions for the weights of the QCS. As long as the number 
of variables (Nqcs) are no fewer than the number of 
equations, solutions always exist. We may also augment 
the system of equations by /q + 1 equations to ensure the 
correct normalization of each codeword. Hence, for N 
which satisfy 

t 

l + lo + loY. ^(^' ™) ^ ^Qcs < V{N/{t -f 1), m) , (5.3) 

codes with m modes correcting t errors exist. 
For example, when m = 2, Eq.(p^) becomes 



lo{t + 2){t+l) 



+ k + l< 



N 
t + l 



1 . 



(5.4) 



which gives a scaling law N k, Note that this up- 

per bound is generally much larger than necessary, as can 
be seen in the examples for i = 3 or f = 4. Much more 
efficient codes may be obtained, because the QCS may be 
chosen to give redundant constraint equations. This may 
be accomplished either systematically (next section), or 
by numerical search (Section VHl). 



VI. CONSTRUCTION ALGORITHM FOR T < 2 
BALANCED CODES 



So far, we have established criteria for and proved the 
existence of bosonic codes for amplitude damping. We 
now present an explicit procedure which obtain a class of 
balanced codes to correct for t = 1 and t = 2 errors. 
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To correct t = 1 error, consider ordered m-tuples 
{xi,X2, ■ ■ ■ ,Xm) such that xi+X2 + . ■ ■ + Xm = n. We wiU 
use the same symbol Q{n, m) for the space of all such m- 
tuples as well as the space of all QCS {|a:ia;2 . . - Xm)}- 
We define an operation, TZ on Q(n, m) which takes 
(xi, X2, . . . , Xm) to (x2, • . • , Xm, Xi), i.e., the symbols are 
cyclically permuted. Define the order of an element to 
be the size of its orbit under TZ. It follows that the order 
p must divide m, and let m ^ pq. An element of order p 
looks like (xi, . . . , Xp, xi, . . . , Xp, . . . , xi, . . . , Xp) with the 
string {xi, . . . , Xp) repeated q times. The orbit looks like: 



(xi , X2 1 ... 1 Xp : 5 ; i ■ ■ • i Xp) 

{X2,X3,...,XI, ,X2,X3, . . . ,Xi) 



{Xp , , . . . , Xp—i , , Xp , Xi , . . . , Xp—i ) . 



(6.1) 
(6.2) 
(6.3) 
(6.4) 



Wc form states by taking equal-weight superposition 
of QCS in each orbit: 

|c) = -^(1x1X2 - --Xp X1X2 - --Xp) (6.5) 

+ \X2X3---Xl X2X3---Xl) (6.6) 

+••• (6.7) 
+ \xpXi---Xp-i XpXi ■ ■ ■ Xp-i)) . (6.8) 

States formed by distinct orbits are orthogonal, as the 
orbits partition Q{n, m). Furthermore, we multiply each 
number in the QCS by d. The minimal separation of 
distinct QCS will be at least c?, since distances come as 
multiples of d only. Hence, all the codewords will remain 
orthogonal after errors oi t < d occur, codewords are 
now in the form: 



-^{\dxidx2 ■ ■ ■ dxp ■ ■ ■ dxidx2 ■ ■ ■ 



\dx2dx3 ■ ■ ■ dxi ■ ■ ■ dx2dx3 ■ ■ ■ dxi 
\dxpdxi ■ • • dxp-i 



(6.9) 
(6.10) 
(6.11) 

■ dxpdxi ■ ■ ■ dxp-ij) . (6.12) 



For the non-deformation criteria, the row sum is nd — 
N by construction. The column sum divided by the nor- 
malization factor squared is (by Lemma 1.1) 



dxi 



P 



dn 
m 



(6.13) 



in any codeword, independent of the order of the con- 
stituent QCS. Codes in examples (l)-(3) in Section VIII 
are constructed in this way. 

To correct for t — 2 errors, the t = 1 criteria have to be 
satisfied as well. We take a subset of the t — 1 codewords 
which will survive the extra non-deformation criteria for 
t = 2. We also replace d = 2 by d > 3. For to > 2, pairs 
of codewords in the form 



-^{\dxidx2- 


• dx^a) 


(6.14) 


+ \dx2dx3 ■ 


■ ■ dxi) 


(6.15) 


1 

-j_ ... 




(P. 1 (\\ 


\ \ m ""^ 1 


■ ■ - dx 


l'6 1 T\ 




dx2dx\) 


f6 18) 


+ \dXrn-l 


■ ■ dxidxm) 


(6.19) 


+ ■■■ 




(6.20) 


+ \dxidxm 


■ ■ ■dx2)) , 


(6.21) 



will always satisfy the non-deformation criteria for t = 2 
(proof omitted). Examples (4) and (5) in Section VIII 
are constructed in this way. This encodes only one qubit; 
we are still looking for t = 2 codes which can encode more 
qubits. 

For t > 3, we performed a numerical search for spe- 
cial QCS in which the system of linear equations for the 
weights is linearly dependent. In the best case, the num- 
ber of linear equations to be solved can be much reduced. 
Therefore we can find codewords involving fewer QCS, 
fewer number of modes and fewer number of quanta. Al- 
though encoding is certainly possible with a much smaller 
Hilbert space, we have not found a systematic way to 
generate such QCS. Cod es correcting t < A errors are 
exhibited in Section VIII 



VII. RATES AND FIDELITIES 

The performance of these bosonic quantum codes can 
be characterized by their rate - number of qubits com- 
municated per qubit transmitted, and by their fidelity 
- the worst-case qubit degradation after decoding and 
correction. We discuss these two measures here. 

The rate r is given by the ratio of the number of en- 
coded qubits to the maximum number of qubits that can 
be accommodated in our Hilbert space: 



k 



TOlog2(n(i + 1) 



(7.1) 



where 2'^=number of codewords, and {nd + 1)™ is the 
size of the Hilbert space in our code. The exact number 
of possible codewords depends on the choice of N and 
TO (we have worked out a counting scheme, but omit the 
details here). However, the majority of the QCS have 
order to. Hence, to a good approximation, the number 
of codewords obtained is: 



2fe ^ V{n,fn) 



(7.2) 



For small n, codewords involving fewer than m QCS al- 
low slightly more qubits to be encoded compared with 
(7.2). This small gain can be important in applications 
such as key distribution in quantum cryptography. 
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We now turn to the code fidelity T , which we desire 
to know as a function of the parameters TV, m and t. It 
is defined as 



J- — min 



(7.3) 



Let the input state be expressed as a sum of codewords, 
IV'in) = "^'Ic/)- Then using the orthogonaUty and non- 
deformation conditions, we find that 



(V'»i4A.iV'm) = i(Qi4AiQ)i, 



(7.4) 



with |q) any one of the codewords in the right hand side. 
Now, if we write each codeword as 



Ml|"-ll"-12 ■ ■ -"Iri 
7i2|f^2lf^22 • ■ ■ "•2ri 



then, for k = {ki,k2, ■ ■ ■ ,km) G /C(s): 



(7.5) 
(7.6) 
(7.7) 
(7.8) 



{ci\AlA~^\ci)\ = {\~^f-^^^ya,C{n,„k,), {7.9) 



1=1 



and using the following relation for binomial cocfhcicnts: 

N, 



C{N,s)= ^ //,C (n,i , fei )C (n^2 , fc2 ) 



we find that the fidelity is 
t 

s=l 

= 1 - C{N, t + l)7*+i + 0(7*+2) 



(7.10) 

(7.11) 
(7.12) 



This expression holds for balanced codes as well as un- 
balanced codes. The amazing feature is that given a code 
which satisfies the orthogonality and non-deformation 
constraints, J- is independent of m; it is determined only 
by N and t. 

One should note that although codes can be con- 
structed to correct an arbitrary number of photon loss, 
the more errors one wishes to correct, the more photons 
are required. On the other hand, having larger photon 
number states means a higher probability for the system 
as a whole to suffer loss of quanta. These two effects 
compete against each other to give an upper bound on 
the fidelity, which we can estimate as follows. Let N the 
required photon number, and t be the total number of 
errors to be corrected. As previously discussed, due to 
the constraint equations which hold for error correction 



to be possible, the two parameters can be reduced to one 
degree of freedom, in terms of which we may estimate the 
optimal achievable fidelity. In terms of t, the optimum 
fidelity for fixed 7 is obtained by setting 



dt 



ln(l - JT) = . 



From Eq.(7.12), this gives to first order in 7 



1 dC dN 1 dC 



(7.13) 



(7.14) 



where C is a short hand notation for C{N,t -1-1). Us- 
ing the Stirling approximation for the factorials in C, we 
obtain 



In 



N \ dN /N -t-1 

N-t-i) ir^ ^ \ t+i 



ln7 = 0. 

(7.15) 



In general, N is much larger than i, which allows further 

ln7 = 0. (7.16) 



simplification of Eq.(7.1f:) 
dN t , f t 



dt N ^^^\N 



Generally, N will asymptotically follow a power scaling 
law in i, so we may approximate N « flat"' where the 
prefactor / and exponent a are approximately constant. 
Solving for the optimum t we find that 



opt 



1/ic 



(7.17) 



Plugging back into the Eq.( [7.12 ) gives us an estimate 
for the optimal achievable fidelity. We thus have a loose 
bound on the fidelity of a bosonic code with arbitrary 
QCS. There is no theoretical bound on the number of 
correctable errors. 



VIII. EXPLICIT CODES 

Some explicit codes resulting from our work are pre- 
sented here. States are given unnormalized when the nor- 
malization factor is common for all codewords. Codes are 
specified as [[iV, m, 2*^, c?]], where N is the total number 
of excitations in the QCS, m is the number of modes for 
each QCS, 2^^ is the number of codewords and d is the 
minimal distance between codewords. The fidelity of all 
the codes are given hy T ^ 1 — C{N, t + 1)7*+-'^. 

Example (1) - [[4,2,2,2]], n = 2, t = 1, fidelity 
T^l- 672: 



|0l)-^[|40) + |04)] 
= |22) 



(8.1) 
(8.2) 
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Example (2) - [[12,3,10,2]], n = 6, i = 1, fidelity 
!F !v 1 — 667^, labels given in hexadecimal (c = 12, 
a= 10): 



Cl) = 


^[|00c) + 


cOO) + |OcO)] 


(8.3) 


C2) = 


^[|02a) + 


|a02) + |2a0)] 


(8.4) 


C3> = 


^[|048) + 


|804) + |480)] 


(8.5) 


C4) = 


^[|066) + 


|606) + |660)] 


(8.6) 


C5) = 


i=[|084) + 


|408) + |840)] 


(8.7) 


ce) = 


;^[|0a2) + 


|20a) + |a20)] 


(8.8) 


C7) = 


;^[|228) + 


|822) + |282)] 


(8.9) 


Cg) = 


;^[|246) + 


|624) + |462)] 


(8.10) 


cg) = 

Ca) = 


^[|264) + 
|444) 


|642) + |264)] 


(8.11) 
(8.12) 



Example (3) - [[6,3,4,2]], n = 3 = 0+ l + 2, f=l, 
fidelity ^ w 1 - 157^: 

|ci) = |600) + |060) + |006) (8.13) 

|C2) = |420) + |204) + |042) (8.14) 

|C3) = |240) + |402) + [024) (8.15) 

\ci) = |222) . (8.16) 

Example (4) - [[9,3,2,3]], n = 3, t = 2, fidelity 
« 1 — 847^: Note this code differs from the previ- 
ous one from having d = 3 instead of d = 2. We take 
only |c2) and jcs) as codewords. 



|0l) = |306) + |063) + |630) 
|U) = 1 036) + 1 360) + 1 603) 



(8.17) 
(8.18) 



Exapmle (5) - [[6,4,2,2]], n=6=0+l+2+3, fidelity 
~ 1 — 157^: The minimal distance between QCS is 
d = 2. However, the QCS are not generated by multiply- 
ing each number by = 2. 

\0l} = |0321) + |1032) + |2103) + |3210) (8.19) 
|U) = |0123) -I- 11230) + ]2301) -|- 13012) . (8.20) 

Example (6) - [[7, 2, 2, 2]], fidelity « l - 217^: The 
codewords are not formed by cyclic permutations of the 
QCS. Note that column one and two have different col- 
umn sums. 



I0i) = l70) + 116) 
1U) = 152) + 134) 



(8.21) 
(8.22) 



Example (7) - [[9,2,2,3]], fidelity T 1 - 84-/^: Un- 
balanced code that will tolerate t = 2 errors. Note that 
one codeword is formed from the other by reversing the 
order of the modes, (this symmetry between the two 
modes is a sufficient condition for balanced codes with 
t = 2, m > 3). 

V3, 



|0L) = i|90) + ^l36) 
|U) = ^|09) + ^|63) 



(8.23) 
(8.24) 



Example (8) - [[9, 3, 2, 3]], fidelity J" « 1-847^: Unbal- 
anced code that will tolerate t ^ 2 errors, showing that 
the symmetry is not a necessary condition for correcting 
t = 2 errors. 



|0l) 



1306) + [360)] 



|li) = ^[V^|333) + ^]009) + ]090)] 



(8.25) 
(8.26) 



Example (9) [[16,2,2,4]], fidelity T ^ 1 - 18207^: 
Unbalanced code that will tolerate t = 3 errors. Labels 
are given in base 17. c and g denote 12 and 16 respec- 
tively. 



\0L) = ^[\0g) + \g0)+VQ\88)] 
lU) = i=[|4c) + |c4)] 



(8.27) 
(8.28) 



Example (10) - [[20,3,2,4]], fidelity J" ^ 1 - 48457^: 
Another unbalanced code that will tolerate t = 3 errors. 
Labels are given in base 21. c, g and k denote 12, 16, 
and 20 respectively. 



i[104.g) +2|40.g) -h2V5lOA;0)] 



|0l) 



= 4e[V2|44c)+V31488)] 
V5 



(8.29) 
(8.30) 



Example (11) - [[50,2,2,5]], fidelity T 1 - 
21187607^: Note the rapid growth in the numerical factor 
in the second term. To correct for large number of errors, 
we need to encode a qubit in a large Hilbcrt space, but 
emission probabilities are large for high number states. 
This puts a limit of performance in our codes. The ac- 
tual code involves numbers five times the numbers shown 
below, a denotes 10. 



|U) 



+ 



I0a) + 

182)- 
119) 
155) 



146) 
191) 
128) 
173) + 



laO). 



(8.31) 

(8.32) 
(8.33) 
(8.34) 
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IX. PHYSICAL IMPLEMENTATION 

Encoding and decoding of the codes we have described 
here can be performed in principle using n-photon eigen- 
states, beamsphtters, phase shifters, and Kerr nonhnear 
optical media. We demonstrate for example how states 
for our simplest code may be constructed. We then dis- 
cuss how decoding and correction may be performed. 

Shown in Fig. ^is a quantum circuit which can be used 
to encode a qubit using the |22), |04) + |40) code. Let 
us see how this circuit works by tracing the state at the 
five indicated points. The initial state is \ipo) = |0122). 
Two 50/50 beamsplitters act on the first and second two 
modes of this state to give 



IV'i) 



|10) + |01) 



V2 



%/6(|04) + |40)) - 2|22) 



(9.1) 



Next, a nonlinear optical Kerr medium is used to per- 
form a cross phase modulation between the two middle 
modes. This serves to "label" the |22) state, giving: 



IV'2 



1 



4\/2 



10) + |01)j [V6(|04) + |40)) 
|10)-|01)) [2|22)] 



(9.2) 



A final beamsplitter in the first two modes now serves to 
turn the phase modulation into a detectable amplitude 
difference, 



1 



2V2 



101) 



\/3sin6'(|04) + |40)) + \/2cos0|22) 



|10) [V3cos6l(|04) + |40)) - \/2sin6l|22)j 1 (9.3) 



such that if the first two modes are measured to be |01) 
(otherwise, the state is discarded) then we have the out- 
put 



re-encoding the state for further transmission. In this 
system, we have a non-binary state. If it can be turned 
into a binary state, then the entire correction procedure 
may be performed using standard techniques of quantum 
computation, with the usual binary quantum logic gates 
p^ . This is possible as follows. 

Consider the circuit shown in Fig. ^ The structure is 
identical with the well-known quantum-optical Fredkin 
gate ||l^,|l3 , but let us think of it here in a different way. 
The lower pair of wires may be considered to be a sin- 
gle "dual-rail" qubit ||l|l, with logical states |0l) = |10) 
and lU) = |01). Wfth an input state \ipo) = |0n)|0L) 
and a Kerr medium with x — '"'/^i then the output state 
will be I'l/'a) = |On)|lL), and in this manner the circuit 
may be thought of as a kind of controlled-not gate which 
distinguishes between a control state of |0) and \n). 

We now use this bosonic controlled-not gate to con- 
struct the circuit of Fig. ^, which is based on the fact 
that different values of x allow us to distinguish differ- 
ent values of n. Furthermore, since the decomposition of 
the number n into sums of powers of two is unique, it 
is convenient to take x to be binary fractions of tt. In 
this manner, the first dual-rail qubit becomes the least- 
significant bit of n, and so-on. If the qubits are measured, 
this circuit would then function cquivalently to a perfect 
photon number detector. 

However, it is much more useful in that calculations 
may be performed based on the binary representation of 
n to determine if any error occurred and to calculate the 
error syndrome. The circuit is then applied in reverse to 
undo the entanglement with the bosonic state, and the 
appropriate correction procedure is applied to fix the de- 
, tected error. This is possible since generalized measure- 
ments may be performed on the qubit states to determine 
the error syndrome without destroying the superposition 
state of the original qubit encoded in the bosonic state. 



IV'4) = cos6l'|22) + sin 61' 



,|04 ) + |40) 
V2 



where 



0' = tan-i 



tan( 



(9.4) 



(9.5) 



is the new effective angle. IV'4) is the desired encoded 
state. After transmission of IV'4) through a lossy com- 
munication link, the final state can then be measured 
immediately (using, for example, photon number coun- 
ters) and the transmitted qubit collapsed. This would 
be the standard procedure for point-to-point quantum 
cryptography. 

Alternatively, the qubit may be relayed by perform- 
ing an error correction step. This involves calculation of 
the error syndrome, correcting any detected error, then 



X. CONCLUSION 

Our treatment of amplitude damping errors is some- 
what unusual from the standpoint of most quantum error 
correction theories, which deal with a bit fiip and phase 
flip picture of errors. The relationship can be understood 
by expressing the ^0 and Ai operators as coherent su- 
perpositions of such errors; from Eqs.(2.13) and (2.14), 



^0 



1 



2 
2 



(1 + + (1 - ^^-l)<Jz 



(10.1) 
(10.2) 



With probabilities up to 0(7), a binary code with m bits 
will either project Af™'\ci) onto a state with no errors, or 
project Af™'~'^Ai\ci) onto a state with only one bit fiip 
error, resulting from a combination of J®"*-! and one of 
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the Pauli operators Ux, <Jy Hence, a binary code correct- 
ing for any one bit error will indeed correct all amplitude 
damping errors up to losing one photon, although not to 
all orders. One reason we have studied bosonic codes is 
to exploit the possibilities for achieving higher efficien- 
cies or easier physical implementation, though the study 
is theoretically interesting on its own. 

It is important to realize that amplitude damping er- 
rors are not independent bit errors, since the decay factor 
of each QCS depends on the total number of excitations 
in it. This fact is also pointed out by Plenio et. al. 
Any code correcting a general one bit error can do so 
in the presence of possible amplitude damping, and is 
capable of correcting any one bit error in the projected 
space in which no loss occurs as well as correcting a single 
photon loss in the absence of any other errors. Two or 
more amplitude damping errors projected into the non- 
identity space are viewed as separate errors. The general 
relation between binary codes and amplitude damping is 
still under investigation. 

Nevertheless, some interesting comparisons may be 
made. Rates from our bosonic codes contrast with those 
achievable by the usual binary codes. For the code of Ex- 
ample (1), = 4, m = 2 and k — 1, so the rate is found 
to be r = 0.22. This is slightly better than r = 0.20 
for the five bit (t, fc) = (1,1) binary perfect code [ p^ , 
and much better than the eight bit (1,1) code of Plenio 
et. al. ^ which corrects errors in the presence of A^'s 
in /C(0). Similarly, for the code of example (2), d = 2, 
n — 6, m — 3 and 2*^ = 10 codewords may be found, 
giving a rate r — 0.2994. In comparison, a naive eval- 
uation of the quantum Hamming bound [ pO[ for binary 
codes gives a possible rate of 0.41. Non-deformation con- 
straints are more restrictive on bosonic error correction 
codes than on binary codes, but the bosonic states admit 
coding schemes which are impossible with binary codes. 
There is no conclusive statement on comparing the gen- 
eral efficiencies of the two different type of codes, but the 
examples we have discovered indicate the existence of a 
rich variety of bosonic codes which may be useful in the 
future. 

The code fidelities may also be compared. Our 
[[4, 2, 2, 2]] code achieves « 1 - 67^. In comparison, 
from an explicit evaluation of the effect of amplitude 
damping on all qubits, we have found that the five-bit 
(1,1) binary code achieves fidelity ~ 1 — 1.757^, while the 
eight-bit (1,1) code achieves only « 1 — 67^! This agree- 
ment with the bosonic code is not accidental; it stems 
from the use of the same total excitation number. How- 
ever, it is worthwhile to point out that despite the effort 
to balance the codewords, the five-bit code still has bet- 
ter performance on average, due to the small number of 
excitation involved in the system. 

In conclusion, we have given general criteria for an 
error correction code which encodes qubits in bosonic 
states. This is a generalization of the binary error cor- 



rection codes. Motivated by the dominant decoherence 
process (amplitude damping) of system such as photons 
transmitted through optical fibers we classify our errors 
according to the number of excitations lost, instead of 
the more common classification of the number of bits or 
modes corrupted. We have shown, in one case, specializa- 
tion to correction amplitude damping does improve the 
rate, number of encoded qubits to number of required 
qubits. However, bosonic codes under amplitude damp- 
ing suffer constraints involving deformation of the Hilbert 
space not shared by the binary codes, rendering the effi- 
ciencies lower in the bosonic case when many qubits are 
encoded. 

It is too early to conclude on the relative performance 
of the binary codes and the bosonic codes. Further study 
will aim at improving the efficiencies, perhaps by using 
relative phases to maintain orthogonality of the code- 
words instead of by using distinct QCS, so that the QCS 
can occur in more than one codeword. Another possi- 
bility is to encode qubits using QCS of different total 
number of excitations. We hope that further study of 
bosonic codes will lead to their practical utilization in 
addition to the current theoretical interest. 
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APPENDIX A: CRITERIA FOR 
NON-DEFORMATION OF HILBERT SPACE 



Consider /C(2) = {0 • • • 02, • • • 20, . . . , 20 • • • 0} IJ 
{0 • • • Oil, • • • 110, • • • , 110 • • • 0} for t = 2 errors. For 
instance, 



^0---02|q) 

^ m 



(Al) 
(A2) 



C(n,^, 2)7^(1-7)^-2 , 
> , \ [nil . . .n™ - 2) , (A3) 



where C(n, m) is the usual binomial coefficient. The norm 
square of this state is 

72(1-7)^-2 " 

(Ci|^0...02^0--02|Q) = 2^nim{nim - 1) • 



2Ni 



i=l 



The term linear in 7 is independent of I (codeword in- 
dependent) by the criteria for t = I; hence, it follows 
that: 



.\^r,2 



2 — 1 



(A5) 



has to be independent of I if the non-deformation criteria 
is to be satisfied. Other A^. with fc = • • • 02, . . . , 20 • • • 
impose the above requirement on other columns. 
Similarly, fc = • • • 11 changes the codeword to: 



/ n,^_in,„72(l - 7)^ 2 

/ ^ \ ^7 \nil . . . Uim-l - lUirn -1) 



Ni 



which has norm square 



7^(1-7)^-^^ 



(A6) 



(A7) 



i=l 



Eq.(4.2) requires the following to be independent of I: 

^ m 



A similar result is obtained for other k with I's at any 
two modes ji and j2- When we allow ji — j2, we include 
the previous result for two photon loss at one mode. This 
proves Lemma 2. 

For arbitrary i, we get equations involving various bi- 
nomial coefficients. Using requirements involving prod- 
ucts of fewer than t riij^s, we can replace the products of 
the binomial coefficients to products involving exactly t 
n.ij . By mathematical induction, the result for arbitrary t 
(A4) is then obtained. This completes the proof of Theorem 1. 
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I out) 



^2) 



IV'3) IV'4) 



FIG. 1. Quantum circuit to encode a qubit using the code of Example (1). As in jl^, signals travel from left to right, 
wires represent optical modes, diamonds represent beamsplitters, meters ideal photon counters, and the Kerr device is an ideal 
nonlinear optical medium which effects cross- phase modulation of angle 7r/2 between single photons. 



|0). 
|n). 



ID 



90 



kerr 



|0) 




90° 





1^0) IV"!) 



1^-2) 



FIG. 2. (left) Optical quantum logic gate used as a building block in the decoding procedure; example input states are shown 
on the left, x is the cross-phase modulation strength of the Kerr medium, which performs the transformation exp(ixa^a6^fe) 
on two modes a and h . (right) Shorthand notation for this circuit. 



n kI2 n/4 




-nil 



I 



|0> 



-7l/4 



-Till 



BitO 
Bit1 
Bit 2 



FIG. 3. Quantum circuit used to decode a bosonic state into qubits. Each thick wire represents a pair of bosonic modes. 
The top wire carries the bosonic state, and the remaining wires are prepared as dual-rail quantum bits, which carry just one 
photon in each pair. 
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